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PROCEDE DE MESURE DE PARAMETRES PHYSIQUES AU MOINS 
UNE PHASE AUX DIMENSIONS MICROMETRIQUES OU 
NANOMETRIQUES DANS UN SYSTEME COMPOSITE. 

5 La pr^sente invention a pour objet un proc£d6 de d^tennination d'au 

moins un param&tre m^canique d*au moins un mat^riau dans un systdme composite 
comprenant au moins deux phases distinctes. Le proc^d6 s'applique plus 
particuli&rement au cas oil la phase k caract^riser est stractur6e k T^chelle 
microm6trique ou nanom^trique, et au moins un autre mat^riau constitue xm substrat 
10 ou une matrice. Par exemple, le proc6d6 s*applique k la caract6risation d'une couche 
mince d^pos6e sur un substrat ou dMnclusions, filaments ou fibres dans une matrice. 

L* expression « param^tres mteaniques » est utilis6e ici pour indiquer 
aussi bien des propri6t6s m^caniques et thermom^caniques, telles que le module 
d'Young, le coefficient de Poisson, le coefficient de dilatation themfiique etc., que des 
1 5 6tats de ddformation et/ou de contrainte. 

La connaissance des propri6t6s m6caniques d'un systfeme ou 
dispositif physique permet d'en optimiser le fonctionnement. Tout mat6riau quel qu'il 
soit est soumis h des contraintes extemes et il faut pouvoir connaitre sa resistance a de 
telles contraintes. II est done important de connsdtre au mieux les propri^t^s 61astiques 
20 d'un syst&me et en particulier celles d'une ou plusiews couches dispos^es a la 
surface d'un substrat. Les propri^t^s de cette couche (ou de ces couches consid6r6es 
comme un sous-systfeme) diffi^rent en effet notablement des propri6t6s des couches 
prises individuellement. 

Les propri^t^s dlastiques des mat^riaux interviennent dans de 
25 nombreux domaines d'applications : rev6tement de pieces m^caniques, deformation de 
structure, ... EUes jouent un role de plus en plus important dans la realisation de 
circuits eiectroniques : la taille des dispositifs diminuant, les contraintes generics a 
rinterface des diffSrentes parties du dispositif sont d'autant plus importantes. 

Les m6thodes d6ji d^veloppees pour mesurer les constantes 
30 d'61asticit6 et les coefficients de dilatation sont g6n6ralement r^alisees k une 6chelle 
macroscopique. 

Les m^thodes traditionnelles sont : 

A) deflexion, deformation d'6prouvettes macroscopiques ou de taille 
plus reduite : voir par exemple : Measurement of elastic modulus, Poisson ratio, and 
35 coefficient of thermal expansion of on-wafer submicron films. Zhao, Jie-Hua; Ryan, 
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Todd; Ho, Paul S.; McKerrow, Andrew J.; Shih, Wei-Yan. PRC/MER2.206, 
Journal of Applied Physics (1999), 85(9), 6421-6424, 

ou bien: Elastic modulus measurements by three methods on a 
2-D laminated continuous fiber-reinforced ceramic matrix composite. Gonczy, S. T.; 
S Jenkins, M. O Ceramic Transactions (2000), 103(Advances in Ceramic-Matrix 
Composites V), 541-547, 

ou bien encore: Analysis of residual stress in cubic boron nitride 
thin films using micromachined cantilever beams. Cardinale, G. P.; Howitt, D. G.; 
McCarty, K. P.; Medlin, D. L.; Mirkarimi, P. B.; Moody, N. R. Diamond and Related 
10 Materials (1996), 5(11), 1295-1302 

B) Ondes ^lectroacoustiques : par exemple Hardness and Young's 
modulus of high-quality cubic boron nitride films grown by chemical vapor 
deposition. Jiang, X.; Philip, J.; Zhang, W. J.; Hess, P.; Matsumoto, S Journal of 
Applied Physics (2003), 93(3), 1515-1519, ou bien Thin-film elastic-property 

15 measurements with laser-ultrasonic SAW spectrometry. Hurley, D. C. ; Tewary, V. 
K.; Richards, A. J Thin Solid Films (2001), 398-399 326-330 

C) nanoindentation par exemple : Determination of mechanical film 
properties of a bilayer system due to elastic indentation measurements with a spherical 
indenter. 

20 Chudoba, T.; Schwarzer, N.; Richter, F.; Beck, U. Thm Solid Fihns 

(2000), 377-378 366-372, ou bien. Hardness and elastic modulus measurements of 
AIN and TiN sub-micron thin films using the continuous stiffness measurement 
technique with FEM analysis. Rawdanowicz, T. A.; Sankar, J.; Narayan, J.; Godbole, 
V. Materials Research Society Symposium Proceedings (2000), 594(Thin Films— 

25 Stresses and Mechanical Properties Vin), 507-512. 

II existe 6galement des methodes plus r^centes pour mesurer les 
propri^t6s m^caniques de petits objets, par exemple la microscopic k force atomique 
('AFM*) d^crite par exemple dans : 

Measurements of elastic properties of ultra-thin diamond-like carbon 

30 coatings using atomic force acoustic microscopy. Amelio, S.; Goldade, A. V.; Rabe, 
U.; Scherer, V.; Bhushan, B.; Arnold, W Thm Solid Fihns (2001), 392(1), 75-84. 

L'invention a pour but de mesurer un ou plusieurs paramfetres 
m^caniques, notamment les constantes 61astiques et les coefficients de dilatation 
thermique d*un mat6riau et les contraintes que ce mat^riau g^nere lorsqu'il est associ6 

35 a d*autres. Sa specificity est de pouvoir mesurer de tels paramfetres dans des conditions 
proches de leur condition d'applications : couches minces ou ^paisses, couches 
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inhomog&nes, ^couches* discontinues ou pr£cipit6s/inclusions, fibres ou filaments, 
boltes (terme plus sp^cifique aux mat^riaux semi-conducteurs oil les charges sont 
localis^es dans ces boftes constitutes d*un second mattriau different du substrat). Les 
connaissances actuelles ne permettent pas d'affmner que les paramitres mesurts 

5 macroscopiquement peuvent £tre utilises dans des structures de petites dimensions, 
c*est k dire micromttriques et nanomttriques. 

L'idte de base de Tinvention consiste k utiliser un ph6nom&ne bien 
connu, mais qui est gtntralement considtrt comme un inconvenient que Ton cherche 
k tviter ou k ntgliger : dans un tchantillon prtsentant au moins une dimension 

10 sufBsanunent faible, par exemple k r^chelle micromttrique ou sub-micromdtrique, du 
fait de la proximity de la surface libre, les contraintes accumul6es dans le dispositif 
initial sont relaxtes (voir les documents prtcitts de Perovic, D., Weatherly G. and D. 
Houghton Phil. Mag. 64 (1991)). Selon Tinvention, il n'est pas n^cessaire d'effectuer 
des mesures directes sur le mattriau k caracttriser (ce qui n'est pas toujours possible 

15 ou ais6) : la mesure des paramfetres de deformation du substrat ou matrice, dont 
Torigine est le relachement des contraintes dans r^chantillon, permet de remonter aux 
param^tres de la phase structurte k rtchelle micrometrique ou nanomttrique. Des 
mesures tventuelles effectutes sur cette phase, lorsqu'elles sont techniquement 
possibles, apportent un complement d*information mais ne sont pas, en general, 

20 indispensables. 

Par « paramdtre de deformation » on entend aussi bien les 
deformations pures, couramment indiqutes par le symbole 8 et qui, du point de vue 
microscopique, correspondent k des variations des param^tres cristallins, que les 
rotations locales des axes cristallins, indiquees par le symbole p. 
25 L'invention porte done sur un precede pour la determination d*au 

moins un parametre mecanique d'au moins un materiau dans un systdme composite 
comprenant au moins deux phases distinctes, caracterise en ce quMl comporte : 

(a) la realisation d'au moins \m echantillon comportant une 
premiere partie d'une premiere phase et une deuxi^me partie d'vme deuxifeme phase, la 

30 deuxieme partie etant constituee du materiau k caracteriser, rechantillon ayant au 
moins une dimension suffisamment faible pour permettre un relSchement des 
contraintes dans ledit echantillon ; 

(b) la mesure, sur ledit echantillon d'au moins un parametre de 
deformation d'au moins ladite premiere phase, en correspondance d'une plxiralite de 

35 points situes k des distances dififerentes d'une interface entre lesdites premiere et 
deuxieme phase ; et 
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(c) la d^teimination, k partir au moins dudit param^tre de deformation, 
d'au moins un param&tre m^canique de ladite deuxi&me phase. 

Selon des modes de realisation particuliers : 
Ledit proc^de comporte: 
S i. la realisation d'une pluralite d'echantillons qui diffl&rent entre 

eux pour au moins une propriety gtometrique ; 

ii. la mise en oeuvre de retape b. sur chacun desdits echantillons ; 

et 

iii. I'utilisation dans retape c. des mesures effectu^es siir ladite 
1 0 plurality d'echantillons. 

L'etape b. est r^p^tee k au moins deux temperatures differentes 
pour au moins un desdits echantillons. 

L*etape c. comporte : 

i. la modeiisation du relSchement des contraintes dans ledit 
IS echantillon (L) en utilisant une premiere estimation d'au moins une propriete 

mecanique du materiau de ladite deuxiSme phase (B); 

ii. la comparaison des resultats des mesures de Petape b. avec ceux 
de ladite modeiisation ; 

iii. la modification de ladite estimation d'au moins une propriete 
20 mecanique du materiau de ladite deuxidme phase et la reiteration des sous-etapes i. k 

iii. jusqu'& la minimisation de Tecart entre lesdits resultats des mesures et de la 
modeiisation. 

Dans ce contexte, le terme « modeiisation » couvre aussi bien les 
simulations numeriques que les modeies analytiques approches. Avantageusement, on 
25 peiit utiliser une simulation numerique par elements finis. 

Ledit systeme composite est choisi parmi : un substrat 
presentant sur sa surface une couche continue, un substrat presentant sur sa surface 
des ilots ou des bandes de metallisation, une couche avec une zone incluse dans le 
substrat, un transistor, une couche a Tinterieur d'un substrat, une matrice contenant 
30 des inclusions, des fibres ou des filaments. 

Ledit echantillon comporte au moins une dimension a rechelle 
micrometrique ou nanometrique. 

Ledit echantillon est une lame ayant deux faces sensiblement 
paralleles et disposees sensiblement peipendiculairement par rapport k Tinterface 
35 entre lesdites premiere et deuxidme phase, auquel cas Petape b. est avantageusement 
repetee pour une pluralite de lames d'epaisseurs differentes. 
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Selon un mode de realisation altematif, ledit ^chantillon est une 
lame dispos^e selon \m angle par rapport & Tinterface entre lesdites premiere et 
deuxi&me phase, auquel cas T^tape b. est avantageusement T6p6t&e pour une plurality 
de lames dispos^es selon des angles dif££rents par rapport k Tinterface entre lesdites 

5 premiere et deuxi&me phase. 

Selon un autre mode de realisation altematif ledit ^chantillon est 
une lame en forme de coin, pr^sentant deux faces formant un angle entre elles, auquel 
cas retape b. est avantageusement r6p6t6e pour une plurality de lames pr^sentant deux 
faces formant entre elles des angles difiT^rents. 

10 - Les mesures pr^vues k T^tape b. sont effectu^es par diffraction 

d'un faisceau 61ectronique convergent. 

L*6tape b. comporte Tobservation de lignes de Holz pour au 
moins un plan cristallographique de ladite premifere phase (A) et la determination d'au 
moins un pararndtre parmi : la largeur des dites lignes de Holz, leur position et leur 

IS stmcture interne. 

L'etape b. comporte la determination d'au moins la largeur d'au 
moins certaines desdites lignes de Holz et le calcul, pour chacune d'entre elles, d'une 
rotation maximale pmax Ic long de Taxe du faisceau d' electrons. 

L'etape c. met en oeuvre le trace d'au moins une courbe 
20 representant une dite rotation maximale en fonction de la distance par rapport k 
rinterface entre lesdites premiere et deuxieme phase. 

Selon un mode de realisation prefere, retape c. met egalement en 
oeuvre par simulation le trace de courbes representant la rotation maximale pmax cn 
fonction de la distance par rapport k Tinterface entre lesdites premiere et deuxieme 
25 phase pour des valeurs possibles du module d'Yoimg et/ou du coefficient de Poisson 
du materiau de ladite deuxieme phase, ainsi que la minimisation de recart entre les 
courbes simuiees et les courbes experimentales pour determiner le module d' Young 
et/ou du coefficient de Poisson du materiau de ladite deuxieme phase. 

L'invention sera mieux comprise k la lecture de la description ci- 
30 apres, en liaison avec les dessins dans lesquels : 

- les figures la i If illustrent differentes configurations du dispositif a 
etudier : couche simple pour la figure la, bandes de metallisations ou ilots pour la 
figure lb, couche avec zone incluse dans le substrat pour la figure Ic, transistor pom- 
la figure Id, couche mince k Tinterieur d'un substrat pour la figvire le et fibres dans 
35 une matrice pour la figure If. D'une maniere generale, on indique par A le substrat ou 
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matrice et par B la phase aux dimensions microm^triques ou nanom^triques k 
caract^riser. 

La figure 2a et 2b reprdsentent une vnt en coupe d'un transistor 
ayant un substrat en silicium Si, surmontd d'une couche de NiSi d*dpaisseur 

5 e = 20 nm et d'une couche de Si3N4. Un contact en tungst&ne constitue le contact D de 
drain* La zone analysde (figure 2b) est pointde par la fl&che F. La lame mince 
representee k la figure 2b prdsente une dpaisseur t (que Ton peut faire varier)/Le 
faisceau incident moyen des Electrons pour les cliches "CBED" (direction zo) est pris 
selon Taxe y2 qui fait un angle y avec Taxe yi normal k la lame. Les directions et les 

10 axes cristallographiques correspondant sont illustrds k la figure 2b. CBED est 
Tacronyme de « diffraction d'^lectrons k faisceau convergent »» Convergent Beam 
Electron Diffraction en langue anglaise. 

La figvire 3 a est un montage de photos illustrant cinq des 
digrammes de dif&actions "CBED" choisis parmi une cinquantaine d'images realis^es 

15 effectivement, selon unedroite perpendiculaire k la surface (direction - Z2), k ime 
distance de 155 lun du contact D de drain. Pour chaque diagramme, on indique la 
distance - Z2 (en nanometres) par rapport k la sxuface du substrat (interface avec la 
couche mince) et Tangle A9 mesure, la courbe de la figure 3b repr^sentant revolution 
de Tangle AO calcuie en fonction de -Z2. 

20 - La figure 4a montre un cliche de diffraction « CBED » pris 

suivant \m axe de zone [230]dans le silicium non contraint, k une profondeur 
importante zi = 300 nm (X2 // [3, -2, 0], y2, // [2, 3, 0], Z2 //[OOl]. 

La figure 4b est le cliche 4a sur lequel on a ajuste des lignes de 
Holz positionnees par le logiciel "jems". La position de ces lignes permet de 

25 determiner avec precision Tangle de demi-convergence du faisceau incident et la 
tension d'acceieration (la methode est ici similaire k celle utilisee dans le projet 
"STREAM"). 

La figure 5a est un cliche CBED en axe de zone [230] dans le 
silicium deforme (diagramme « CBED » de type c pris a la profondeur Z2== 140 nm). 
30 Sur la figure 5b, on a supeipose k ce cliche des lignes issues de 

3 simulations avec le logiciel " jems " : 

■ le premier systfeme de lignes represente la simvilation du cliche 
de diffraction d'un cristal de silicium parfait desoriente de +9niax par rapport k Taxe 
X2// [3-20] ; 
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■ le deuxi&me syst^e de lignes repr^sente la simulation du 
cliche de diffraction d*un cristal de silicium parfait d^sorient^ de -Gmax par rapport k 
raxex2//[3-20]; 

■ le troisi&me syst&me de lignes (en pointiUdes) repr^sente la 
5 simulation du cliche de diffraction du cristal parfait non d^sorient^. On note par A6 

Tangle 2emax.. 

- Les figures 6a k 6d illustrent la faisabilit^ de minimiser les 
coefBcients d'61asticit6 du matdriau. 3 paramdtres ont 6X6 minimises manuellement par 
essais/erreurs : la temperature de coherence To (temperature Active oi les deux 
10 mat^riaux seraient cohdrents), le module d' Young E et les coefficients de Poisson v. 

Les carr6s repr^sentent les mesures experimentales de Tangle AO. Les 
courbes avec de nombreux points sent obtenues k partir des simulations par elements 
finis et repr6sentent Tangle 0,98A|3 ou AO 2 pmax- Le critdre de minimisation % 
lvalue la distance entre les angles AO mesures et 0,98Aj5 calculus, k diff^rentes 
15 profondeurs Z2 sous le siliciure et k difKrentes ^paisseurs t de lame. 

La figure 6a : ^paisseur de la lame t = 300 nm, 
valeurs initiales des parametres pris dans la litt^rature : T = 410®C, 
E=150GPa,v = 0,l 

La figure 6b : ^paisseur de la lame t = 300 nm, 
20 valeurs initiales des paramdtres pris dans la litt6rature : To - 410^C, 

E=150GPa,v = 0,l. 

La figure 6c : ^paisseur de la lame t - 320 nm, 
valeiirs finales apr^s une minimisation manuelle 
partielle : To = 430°C, E = 1 15 GPa, v = 0,288. 
25 - La figure 6d : ^paisseur de la lame t = 320 run, 

valeurs finales apr^s une minimisation manuelle 
partielle : To = 430^C, E = 1 1 5 GPa, v - 0,288. 

La figure 7a est un diagranmie de diffraction CBED obtenu 
pour une lame d'6paisse\ir t = 320 nm k la distance -Z2 = 84 nm de la surface. 
30 - La figure 7b est un profil qui a 6X6 r6alis6 sur la bande de Holz 

(5,-3,7) de la figure 7a. L'6cart angulaire A0g est marqu6 par 2 traits verticaux et 
represente la largeur de la bande de Holz. 

La figure 7c est une simulation du profil de la figure 7b qui 
montre qu'il est possible de reproduire T61argissement des lignes de Holz et les 
35 variations d'intensit6 1 - {&) dans les bandes de Holz. Cette simulation a utilise les 
r^sultats du calcul par 616ments finis qui est une illustration du point iv (constantes 
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^lastiques optimis6es et d^placement R(yhZi)). A9g est 6gal k 0»1 1°, ce qui correspond 
' k une rotation d'axe X2 d*angle AO - 0,14^. Dans les calculs pr^liminaires» on n'a pris 
en compte que les largeurs AOg d*axe X2//[320]. Un algorithme de minimisation plus 
perfectionn6 chercherait & reproduire les courbes similaires h 7b dans leur ensemble, 
5 avec leurs oscillations et non pas k reproduire seulement leur largeur AOg k travers 
Tangle AO. 

Les figures SkSd illustrent im deuxi&me exemple experimental, 
dans lequel le substrat A de silicium est surmont6 d'une couche de Si(i.x) Gex, puis 
d'une couche de Si en surface. La figure 8a montre la lame fine d'6paisseur t. Les 

10 figures 8b k 8d sent des diagrammes de diffraction " CBED " respectivement dans la 
direction [230] k rint^rieur de la couche Si(].x) Gex, dans le substrat Si loin de la zone 
d^form^e, et dans le substrat Si pr^s de la zone d^form^e. 

La figure 9 montre un organigramme d'un mode de realisation 
du precede de Tinvention. 

IS La description detaill^e qui suit porte sur le cas particulier oil la 

phase k caract^riser est une couche mince, les ^chantillons utilises ont vme forme de 
lame mince avec des epaisseurs diff^rentes, le paramdtre de deformation consider^ est 
la rotation des axes cristallins dans un plan perpendiculaire k la couche et la technique 
de mesure adoptee est la diffraction d'eiectrons a faisceau convergent (CBED). 

20 Dans ce cas, le proc^de de Tinvention comporte : 

a) la realisation d'une lame L d'epaisseur t suffisanunent faible et 
ayant deux faces sensiblement parall^les et disposees sensiblement 
peipendiculakement k ladite surface de substrat ; 

b) la mesure sur ladite lame d'au moms un paramStre de 
25 deformation du substrat k differentes profondeurs par rapport k la surface; par 

paramdtre de defoimation, on entend egalement parametre de rotation du substrat ; 

c) la determination a partir au moins dudit paramStre de 
deformation/rotation d'au moins un paramfetre mecanique de ladite couche. 

Le procede peut comporter la realisation de plusieurs lames 
30 d'epaisseurs differentes ainsi que la mise en oeuvre de Petape b sxir chacime desdites 
lames. Pour au moins une dite lame, retape b peut 6tre repetee k au moins deux 
temperatures differentes. 

Ladite mesure est avantageusement realisee en generant, pour des 
points du substrat situes a differentes profondeurs, des diagrammes de diffraction d'un 
35 faisceau eiectronique convergent (CBED) d'axe Zq desoriente par rapport k la normale 
a ladite lame, lesdits diagrammes comportant des lignes ou bandes de Holz. La 
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determination c peut alors comporter le lelevd de la largeur des lignes de Holz d*au 
moins certains desdits diagranunes, pour au moins un plan cristallographique du 
substrat. A partir de la largeur de ces lignes de Holz, on peut calculer potir chaque 
diagramme une rotation maximale pmax le long de Taxe du faisceau d' Electrons. Cette 
5 rotation est induite par la couche (ou les couches) disposde sur le substrat et elle en 
caract6rise les propri^t^s* 

On peut ensuite tracer au moins une courbe repr^sentant ladite 
rotation Pmax en fonction de la profondeur k laquelle ont 6t6 obtenus lesdits 
diagranmies. On peut ensuite tracer par simulation au moins une courbe repr6sentant 
10 cette rotation en fonction de la profondeur pour des valeurs possibles des coefficients 
k extraire par exemple le module d' Young et/ou le coefficient de Poisson de la couche 
k analyser dans le cas d*une approximation isotrope et minimiser T^cart entre au 
moins une courbe simul^e et une courbe exp^rimentale correspondante, pour 
determiner ces coefficients. Une technique similaire pourra 6tre utilis^e dans le cas 
15 d'une mod61isation anisotrope avec les coefficients connus de THomme du Metier. On 
prendra pour valeur les param^tres ayant conduit k la courbe simul6e la plus proche de 
la courbe experimentale. 

On pr^sentera ci-aprds les quatre points pr^cit^s de Tinvention. 
point Ci) : La lame mince k g^ometrie contr616e est extraite ou 
20 amincie dans le dispositif. Une lame k face parall^le est preferable mais pas 
indispensable. Un leger angle peut etre present. On a utilise un faisceau d'ions focalise 
TIB', mais des procedes altematifs, classiques en preparation d*echantillons pour la 
microscopic eiectronique peuvent 6tre utilises (amincissement mecanique, clivage,...) 
mais la technique 'FIB' a Tavantage d'etre rapide et de ne pas perturber 
25 mecaniquement le systdme ou dispositif et de contrdler entierement les operations. 

point (ii) : Pour mesurer les deformations, on utilise par exemple un 
faisceau d' electrons convergent (CBED) en mode balayage. Un faisceau eiectronique 
de taille nanometrique est focalise en dififerents points de rechantillon de fagon k 
cartographier les deformations completes de la lame mince (figure 3). En chaque point 
30 un diagramme de diffiaction *CBED' est done obtenu. Ce diagranune comporte de 
nombreuses lignes de Holz (plus de 10 lignes) qui correspondent chacune a un plan 
cristallographique indexe par un vecteur g du reseau reciproque. U est preferable 
d'acquerir les diagrammes 'CBED' dans le substrat (partie a), car la partie b est 
generalement trop fine ou trop defectueuse pour que Ton puisse acquerir de bons 
35 cliches de diffraction dans cette zone b. Mais pour certains systfemes bien cristallises, 
des cliches peuvent Stre obtenus dans la zone a et b. Ceci serait par exemple le cas si 
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la couche a est du silicium et la couche b est compos^e de deux couches fines, une 
couche de Sii-xGcx (composition x faible en Ge, par exemple 10%) et une couche 
de Si. 

Pour une 6p£dsseur de lame donn^e et une temperature d'observation 
S donn6e tout un ensemble de cliches CBED est realise de fa9on k determiner les 
rotations et param^tres cristallins du substrat 

point (iii) : En faisant varier Tdpaisseur de lame (la technique TIB* 
permet de r^aliser cela facilement) et la temperature d'observation dans le microscope, 
on obtient tout un ensemble de donn^es experimentales qui permettront (point iv) de 
10 calculer les constantes d'61asticite, le coefficient de dilatation de la partie B ainsi que 
les contraintes dans le dispositif initial non aminci. Dans le cas oti la partie b est 
composie d'une couche homogdne d'un materiau donne, il n'est pas ndcessaire que les 
mesures soient r^alis^es sur ime m6me lame successivement amincie k difBSrentes 
epaisseurs. Travailler sur une lame unique augmente toutefois la precision et est 
15 indispensable dans le cas ou le syst^me est constitu6 d'un nanosystfeme (transistor par 
exemple) unique. Une seule epaisseur de lame ne permet pas de calculer correctement 
toutes les constantes du materiau* Une seule epaisseur de lame ne donne des 
informations que sur les contraintes du mat^riau (c*est le cas de Vetude de la courbure 
des tranches semiconductrices ou 'wafers' via la formule de Stoney). Les contraintes 
20 sont en partie relfichees par une cowbure du substrat (voir par exemple : Measurement 
of elastic modulus, Poisson ratio, and coefficient of thermal expansion of on*vyrafer 
submicron films. Zhao, Jie-Hua; Ryan, Todd; Ho, Paul S.; McKerrow, Andrew J.; 
Shih, Wei-Yan. Journal of Applied Physics (1999), 85(9), 6421-6424) 

point (iv) : Des simulations utilisant la fheorie eiastique sont ensiiite 
25 effectuees pour reproduire les r^sultats experimentaux. Dans des systdmes complexes 
seuls des calculs par elements finis peuvent 6tre realises. Dans des systemes plus 
simples (partie B = couche mince) une formule analytique peut etre impiementee. Les 
simulations reproduisent les phenom^nes suivants : 

- Dans le dispositif initial, les differentes parties sont contraintes ou 
30 partiellement contraintes - Le fait d'extraire une lame mince du dispositif initial 
relSche les contraintes sous la forme de rotation/changement de paramfetres cristallins 
et c'est Tobservation et la simulation de cette relaxation des contraintes qui nous 
permet de determiner les parametres de la partie B. 

L'invention trouvera de nombreuses applications dans le traitement 
35 de surface de pieces mecaniques, d'optimisation de circiiits eiectroniques (contact 
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m^tallique, couche d'oxyde ...) ou dispositifs oi^ la presence de deux mat^riaux 
diff&rents cr6e nicessairement des contraintes m^caniques. 

Le proc^d6 selon rinvention est original bien quil fasse appel h des 
techniques ou des effets physiques bien connus : 
S - la technique du faisceau d'ions focalis^ ou 'FIB' (Focus Ion 

Beam) pour la preparation d'^chantillon. 

le faisceau convergent (CBED : Convergent Beam Electron 
Diffraction), une technique particuliire k la microscopic 61ectronique. Le projet 
STREAM (voir par exemple la publication 'Software for automation of TEM/CBED 
10 Methodology for strain determination' IST-1 9999-1 0341 STREAM Consortium - 
Deliverable D23) utilise une technique proche (*FIB\ faisceau convergent, 
simulation), mais la technique mesure essentiellement des variations de param^tres 
cristallins alors que le proc6d6 selon Tinvention s'int^resse principalement aux 
rotations locales du r^seau cristallin. De plus le projet STREAM ne cherche pas k 
15 mesurer des constantes 61astiques, mais i mesurer des deformations contraintes dans 
des circmts int6gr6s. 

Le projet STREAM ne mesurait que des variations de paramfetres 
cristallins loin des deux parties A et B du dispositif et n^gligeait la relaxation des 
contraintes dans la lame mince. La detection de rotation du r^eau cristallin selon 
20 rinvention permet d'Stre plus rapide, plus pr6cis et de se rapprocher de I'interface 
entre les parties A et B. 

La puissance du proc6d6 selon I'invention a €t6 montr^e en analysant 
les deformations introduites par ime couche de NiSi dans un circuit int^gr^ 
(6chantillon n^l, figure 2a). Ce systSme est reladvement complexe, car il y a plusieurs 
25 contacts 61ectriques, plusieurs matdriaux presents. Pour Stre trait6 rigoureusement 
(ce qui est tout k fait possible), il faudrait tenir compte de toutes les composantes du 
syst&me. 

Dans une premiere analyse, on a suppose que la partie B pouvait etre 
consid6r6e comme une couche mince de taille lat^rale infinie surmont^e par 
30 Tatmosphere (les couches situ6es au-dessus de la couche NiSi 6taient 
exp^rimentalement d6coll6es, T^paisseur de la couche mince a 6t6 mesurfe 6gale k 
e=20nm (voir figure 2b). 

On d6crira ci-apr6s de fa9on precise la realisation des quatre points 

de la partie 3. 

35 Pour cela, il feut d^finir diflGSrents rep^res g6om6triques. 
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Ro =(xo,yo,2o) d^signe le rqjfere g6om6trique \\6 au microscope. L'axe 
zo est paiallile k Taxe optique du microscope, d6fini comme 6tant la direction 
moyemie selon laquelle les Electrons se propagent avant Tdchantillon. Les images ou 
les cliches de diffraction sont enregistr^s dans le plan (xo,yo)* 
5 On d^fmit 6galement difSSrents repdres Rc=(Xc,yc>2cX Ri=(xi,yi,zi), 

R2Kx29y2>Z2)> H^s k la structure cristalline du substrat de silicium (partie A) : 

xc=[100], y£=[010], 2p=[001] 
. Xj=1/a^[110], y, =1/V2[110], 2,=[001] 
10 X, =1/Vl3[320], y, =l/>/l3[230], Z2=[001] 

Le repfere R2 se d^duit du repfere Ri par une rotation d'axe [001], 
d'angle Y=n.31*'. (voir schema 2b), Les relations entre le repfere du microscope et les 
repferes cristallins dependent de Torientation de T^chantillon dans le microscope. 

15 Point Ci) : La lame mince it faces pratiquement parallfeles k €t6 

r6alis6e avec xm TIB' (voir remarques 1 et 2). La normale aux faces a €t6 choisie trds 
proche de la direction yi, (figure 2b) (mais ici aussi \me g6om6trie diff^rente poxjrrait 
6tre choisie : par exemple y2 ). Dans la premifere s6rie d'exp6rience, I'^paisseur t de la 
lame a 6t^ de t='320nm. Cette ^paisseur a 6t6 mesur^e par une technique relativement 

20 classique de faisceau convergent - voir le livre 'Electron MicroDif&action' de J.CH. 
Spence et J.M Zuo (Plenum Publishing Corporation). 

Remarques : 1) L'^cart de parall61isme a 616 mesur6 par faisceau 
convergent et perte d'^nergie. Ces techniques sont relativement classiques (cf. le livre 
pr^cit^ •Electron MicroDi£Braction' de J.CH. Spence et J.M Zuo pour le faisceau 

25 convergent et le livre TEnergy-Filtering Transmisssion Electron Microscopy' de 
Reimer publi6 chez Springer Verlag pour les pertes d'6nergie). On a done mesur6 
r^paisseur de la lame k diffSrentes distances de la couche sup6rieure. On a trouv^ que 
la lame sup6rieure fait un angle de 1.15** par rapport k la lame inf^rieure. Dans notre 
demonstration de faisabilit^, cet angle a 6t6 n6glig6 et on a consider^ que la lame avait 

30 des faces paralldles. Mais il aurait 6t6 possible d'introduire cet angle dans les calculs 
d'616ments finis. 

2) Le FIB est 6galement capable de produire des lames avec un 
meilleur parall^lisme. 

Point (il) : Les diagrammes 'CBED' ont 6X6 pris dans (ou tres proche 
35 de) la direction d*orientation y2=[230], c'est k dire que y2 est parallfele a zq. Cest la 
direction utilis6e dans le projet STREAM, (mais d'autres directions d'observation sont 
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6galement possibles. L'obtention de diagrammes CBED dans difKrentes directions 
viendrait augmenter le nombre de donn^es exp^rimentales). Typiquement, la taille du 
faisceau ^lectronique a €ti prise ^gale k 0,4 nm. Tangle d'ouverture du faisceau est 
voisine de ISmrad. Lors d'un balayage, de 50 k 100 diagrammes CBED sont pris tous 
S les 4 nm dans une direction perpendiculaire k la surface (direction Z2), en partant de la 
surface, mais seulement une dizaine est retenue dans les calculs. La figure 3 indique la 
position selon I'axe Z2 ot 5 de ces cliches exp^rimentaux ont 6t6 obteniis. 

Plusieurs types de diagrammes CBED peuvent 6tre ddfinis : 

a) Loin de interface entre les parties A et B, le substrat de silicium 
10 est, avec la precision de mesure, consid^r^ comme non d6form6 (parexemple le 

diagramme -Z2=300nm de la figure 3a). II peut servir de r6f6rence pour determiner les 
param^tres exp^rimentaux comme cela est fait dans le projet STREAM. 

b) En se rapprochant de Tinterface, les lignes de Holz restent fmes 
mais commencent k se ddplacer l^g^rement : les paramfetres cristallins du substrat 

15 subissent de trfes 16g6res modifications. C'est Teffet mesur6 et quantifi6 dans le. projet 
STREAM. 

c) Plus pr^s de Tinterface, les lignes de Holz s'61argissent (pea 
exemple le diagramme -za^^nPnm de la figure 3a). Nous les app^lerons bandes de 
Holz. C'est Tefifet qui est mesur^ et quantifi^ ci-apr&s. 
20 d) Encore plus proche de Tinterface, les lignes de Holz deviennent 

trop larges et trop faibles : les diagrammes 'CBED* ne contiennent plus assez de 
details pour gtre quantifies (par exemple le diagranune *Z2=70nm de la figure 3a). 
Point f iii^ : 

Une fois cette s^e de mesures r^alis^e k T^paisseur t=320nm, la 
25 lame a €t6 amincie une nouvelle fois jusqu'i une ^paisseur de t=300nm. Une deuxi&me 
s^rie de mesures a alors ^t^ r^alis^e. 

On a choisi 14 diagrammes pour la premiere 6paisseur et 10 
diagranunes pour la deuxieme (voir figure 5). Ce nombre est suffisant pour reproduire 
les variations de la rotation. 
30 Point (iv^ : 

On a utilis6 le logiciel "jems" (P. Stadelmann, CIME-EPFL 
CH1015-LAUSANNE) pour reproduire la position des lignes de Holz des cliches 
CBED loin de Tinterface (cliches CBED de type a), mais d'autres logiciels ou la 
th6orie pr6sent6e dans le livre pr6cit6 de Spence et Zuo pourraient 6tre utilises. Si Ton 
35 reporte ces lignes de Holz du substrat de silicium — le oistal de r6f6rence - sur les 
cliches CBED de type c, ces lignes se situent au milieu des bandes de Holz (la lame 
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est homog&ne et sym^trique). Les positions d'une des extr^mit^ de ces bandes de 
Holz sont simul€es dans le logiciel "jems" en d^sorientant le crista! de r^fiirence par 
une rotation d'angle -Omax, d'axe X2 // [320] (lignes de la figure Sb)* Les positions de la 
deuxiSme extr^mit^ sont simul^es dans le logiciel "jems" en d^sorientant le cristal de 

5 r^fiSrence par une rotation d'angle 8max> d'axe X2 // [320] (lignes de la figure Sb). Une 
rotation d'angle AO r 20maxpermet done de reproduire avec une bonne approximation, 
les largeurs AOg variables des bandes de Holz. 

Avant Tamincissement par la technique TIB', le substrat de silicium 
est contraint, mais g^n&alement peu d^form^. Aprfes amincissement, les contraintes 

10 engendr^es par Tinterface entre les parties a et b sont rel&ch6es. En chaque point de 
coordonn^e (yi, zi), les 616ments finis donnent done un d^placement R(yi, zi) de 
composantes u et v sxir les axes yi et Zi. Ce ddplacement peut gtre d6compos6 en une 
translation locale, une deformation pure e(yi, zi) et vme rotation locale d'axe xi et 
d'angle P(yi, zi) : 

Les simulations dans le logiciel jems indiquent que ce sont ces 
rotations locales p(yi, zj) le long du faisceau d'dectrons (direction zo) qui sont la 
cause principale de I'^largissement des lignes de Holz en bandes de Holz. Une rotation 
d'axe xi et d'angle p(yi, zi) peut 6tre d^compos^e en 3 rotations d'axe X2, d'angle 
20 9 (yi, Z|), d'axe ya, d'angle 0'(yi, zi) et Z2, d'angle 9"(yi, zi). Au premier ordre, seule la 
premiere de ces 3 rotations est importante et Ton a r6galit6 : 

e(yi, zi)=0,98 P(yi, zi) 

25 En particxilier, on a 6galement montre que Tangle Gmax est au premier 

ordre 6gal k 0,98pmax , oil pmax est la valeur maximale de Tangle p(yi, zi) le long du 
faisceau d'61ectrons, c'est k dire le long de y2, ou ce qui revient au meme k la valeur 
maximale de Tangle P(yi, zi) le long de la direction de yi puisque notre lame est 
p^riodique en Xi : 

30 ^max(^i) =0.98/?„ax(^i) = 0.98 Afax(>ff(j^„2|)) 

Dans cette demonstration de faisabilit6, on utilise cette propri6t6 
approch^e pour optimiser les constantes du mat6riau. Chaque diagranmie CBED pris k 
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une profondeur zz est done caract6ris6 par un seul paramdtre A9(Z2) 20max (zi) 
(figure 3b). 

Les courbes de la figure 6 tracent les valeurs mesur^es 
exp^rimentalement et les valeurs 0,98*^ (avec Ap=2 *Pmtod d6termin£es k partir 
5 des calculs par 61^ments finis. 

Les courbes calcul6es de la figure 6a et 6b utilisent des paramitres du 
mat6riau NiSi pris dans la litt^rature. Pour reproduire au mieux les courbes 
exp6rimentales, on a fait varier les param&tres du mat^riau NiSj dans le calcul des 
courbes calcul6es des figures 6c et 6d. Par essais/erreurs, on obtient les constantes de 
10 la couche NiSi (module d' Young E et coefficient de Poisson v (en se pla9ant dans 
Tapproximation isotrope). 

On donnera maintenant quelques precisions sur I'utilisation des 
elements finis. On a utilise un logiciel du commerce. On s'est plac6 dans 
Tapproximation des deformations planes. 
1 5 Trois 6tats particuliers sont importants. 

6tat 1 : On a suppos6 qu'il existe une temp6rature To=Ti+AT k 
laquelle le substrat de Silicium et la couche de NiSi seraient coh6rents, sans contrainte 
(cette temperature n'est pas forcement accessible, on pourrait la designer conune la 
temperature fictive de coherence*). 
20 etat 2 : La lame mince d'epaisseur t, coherente k la temperature To est 

refiroidie k la temperature experimentale Ti=To-AT. Les contraintes sont en partie 
relaxees k la surface de la lame mince. 

On a utilise les lois classiques de la mecanique des solides, pour 
determiner les deplacements R(yi,zi). En particulier la relation entre contraintes et 
25 deformations est donnee par : 



(l + v)(l«2v) 



1-v V 0 V 
V 1-v 0 
0 0 l-2vj(^ 









yiyi 


E 






a 


(l-2v) 




yi2i> 







AT 



ou AT = To-Ti, E est le module d'Young, v le coefficient de Poisson, 
30 a le coefficient de dilatation thermique. Ces relations s'appliquent dans le substrat de 
Si (EA=156GPa, va=0,277 et aA=2,6e-6 K'^) et dans la couche mince de NiSi (Eb, vb 
ttB et AT ont ete optimises en minimisant la distance % entre les coxirbes 
experimentales et calcuiees (voir figure 6)) en utilisant Thypothese des deformations 
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planes. Les deformations sont calcul^es dans le plan d^fini par la direction des 
Electrons incidents yi // zO et Taxe zi. 

fetat 3 : Avant extraction de la lame, le systime ^tait beaucoup plus 
contraint car les contraintes ne sont pas rel&ch^es par les surfaces, line fois les 

5 param&tres (par exemple le coefBcient d'^lasticit^ et le coefficient de dilatation 
thermique) ^valuds, on peut grfice au present proc^d^, ^valuer les contraintes dans le 
circuit int6gre avant amincissement. Les contraintes dans le transistor sont calculSes 
en utilisant les constantes optimis^es : ce sont les contraintes prises au milieu d*une 
lame de tr^s grande ^paisseur to ou celles d'une lame pSriodique d'^aisseur infinie. 

10 Remarques : 

Rl) Pour le point (iv), dans cette demonstration de faisabilite, on a 
ajuste manuellement les constantes d'61asticite, mais ime optimisation automatique est 
possible. 

R2) Pour le point (iv), nous avons 6galement utilise la relation 
15 approchee : 9max = 0.98 pmax, c'est i dire que Ton n*a pris en compte que la largeur 
totale AOg des bandes de Holz. Un ajustage plus precis mais plus long consisterait k 
reproduire les bandes de Holz dans leur ensemble. Chaque diagramme CBED ne serait 
plus alors caracterise par une seule valeur A9 mais par un ensemble de profils 
identiques & ceux realises k la figure 7b. 
20 Ces profils reprdsentent en fait la fonction l-Ag^(s) o(i Ag(s) est 

Tamplitude de la fonction d'onde difiOractee par les plans d'indice g et s mesure Tecart k 
la loi de Bragg pour ces plans g. Comme explique dans le livre pr^cite de Spence et 
Zuo, s et G sont relies par une formule simple et differentes approximations peuvetit 
6tre realisees pour calculer Ag(s). On a montre que pour les lignes fines de Holz de 
25 diagrammes CBED pris dans la direction y2//[230], une approximation deux ondes ou 
tme approximation cinematique doiment le m§me resultat physique. Dans 
Tapproximation cinematique : 




30 Le parametre k minimiser serait dans ce cas la sonune des distances 

entre des courbes du type de la figure 7b et 7c. L'inter6t de la formule precedente est 
que non seulement les rotations, mais egalement les deformations pures sont prises en 
compte dans le terme R(y2,Z2) (mSme si on a montre que Teffet essentiel dans 
reiargissement des bandes de Holz est celui du aux rotations locales). 
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R3) Les points (i) et (ii) ont €t6 ^galement r^alis^s sur r^hantillon 
d6crit dans la partie 3. 

R4) A difiiSrentes temperatures, le principe de la m^thode reste 
similaire et est relativement classique dans la mesxire des deformations. Introduire des 
5 mesures k diff^rentes temperatures revient k augmenter le nombre de donn^es 
experimentales et le nombre de param&tres k optimiser : la mithodologie est tout k fait 
similaire. 

L'organigranune de la figure 9 donne une vue synthetique du mode 
de realisation du procede de Tinvention qui vient d'Stre decrit en detail. 

10 L'etape El du procede est la decoupe d'au moins une lame 

comprenant une partie du substrat et une partie de la couche disposee k sa surface. 
Cela entraine un rel&chement des contraintes (E2), et done une deformation de ladite 
lame. A Tetape E3, au moins un paramfetre de deformation (de preference une rotation 
p) est mesure en correspondance d'une pluralite de points de la lame, k differentes 

1 5 profondeurs par rapport k I'interface couche-substrat. 

D*\ine maniere independante des etapes El - E3, on effectue ime 
modeiisation, typiquement par elements finis, de la lame contrainte, c'est k dire avant 
sa decoupe. Pour cela, on effectue une premiere estimation des proprietes mecaniques 
de la couche mince, telles que son module d' Young E% son coefScient de Poisson v' 

20 et sa temperature de coherence T'o (etape E4). Ensuite, k retape E5, on modeiise le 
relfichement des contraintes. Cela pent se faire en rempla9ant des conditions au 
contour de deplacement impose (nul) par des conditions de contour libre. A retape E6 
on detemiine, sur la base de cette simulation, les valeurs attendues e', P' du ou des 
parametres de deformation en correspondance des points oil ont ete effectuees les 

25 mesures de retape E3. 

De preference, mais pas obligatoirement, les etapes El - E6 sont 
repetees pour une pluralite de lames d'epaisseurs differentes, ou plus gen&ralement 
presentant une geometric differente, et/ou k des temperatures differentes. 

L*erreur quadratique moyenne entre les valeurs mesurees et attendues 

30 est calcuiee (E7), et minimise par variation des estimations des proprietes mecaniques 
de la couche mmce (E9) et iteration des etapes E4 - E7, jusqu'Jl qu'xrne convergence 
soit atteinte (E8). 

De cette fa9on, on obtient une estimation optimale, au sens de 
Terreur quadratique moyen, des proprietes mecaniques de ladite couche (ElO). 
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Optionnellement, Tinformation ainsi obtenue est utilisie, k T^tape 
El 1» pour calculer un 6tat de d^fonnation et/ou de contrainte de la couche mince et/ou 
du substrat. 

La description ci-dessus a 6t6 essentiellement limit^e au cas d'un 
5 substrat sur lequel a 6X6 d6pos6e une couche mince et k des mesures effectu^es sur des 
6chantillons constitu6s par une lame k face sensiblement parall&les entre elles et 
pexpendiculaires k Tinterface entre le substrat et la couche. Cependant, la m^thode 
s'applique ^galement k des ^chantillons d*une forme diff<^rente, par exemple en forme 
de coin. Ce qui est important est que la g^om^trie de cet dchantillon soit connue avec 
10 precision, pour pouvoir effectuer la mod^lisation, et qu'elle pemiette vai rel&chement 
des contraintes. 

De la mSme fa9on, lorsqu'on effectue des mesures sur une pl\u:alit6 
d'6chantillons pour obtenir une meilleure estimation des propri6t6s m^caniques de la 
couche, on pent faire varier une quelconque propri6t6 m^canique desdits ^chantillons, 

15 et non seulement leurs 6paisseurs. Par exemple, on peut utiliser une s^rie 
d'^chantillons en forme de coin pr6sentant des angles difKrents entre leurs faces et/ou 
une diff^rente inclinaison par rapport k Tinterface couche-substrat. 

Comme on s'est int^ress6 particulidrement k la caract^risation d'une 
couche de surface, les diffiSrentes mesures de paramitres de deformation ont ^t^ 

20 effectu^es k des profondeurs diffi§rentes par rapport k I'interface couche-substrat. Plus 
g^n^ralement, dans le cas d'inclusions, fibres ou filaments, on ne pourra pas parler de 
« profondeur », mais simplement de distance de I'interface entre les deux phases 
consid^rdes. En outre, certains des points de mesure pourront, au moins dans certains 
cas, 6tre situ6s k Tint^rieur de la phase nanom^trique ou microm^trique k caract^riser. 

25 En ce qui conceme la technique de mesure, d'autres m^thodes que la 

CEBD peuvent etre utilis^es, comme par exemple la LACBED (Large Angle 
Convergent Beam Electron Diffraction : diffraction d'un faisceau convergent 
d'61ectrons k grand angle). 

Dans r exemple consid6r6, le seul paramfetre de deformation 

30 consid6re 6tait Tangle de rotation P, d6teraiin6 par la mesure de Telargissement A9 
des lignes de Holz. Le d^placement des lignes de Holz, li6 aux d6fomiations pures, 
aurait 6galement pu etre pris en compte. Plus g^n^ralement, le proc6d6 de Tinvention 
peut comporter Tutilisation d'un ou plusieurs parametres de deformations, determines 
k partir de differentes grandeurs mesur^es directement. 
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REVENDICATIONS 
\. Proc^d^ pour la d6tennination d'au moins un param&tre 
m^canique d*au moins un mat^riau dans un systdme composite comprenant au moins 
deux phases distinctes (A, B), caract^ris^ en ce qu'il comporte : 

S a. la realisation d'au moins un ^chantillon (L) comportant une 

premi&re partie d*une premi&re phase (A) et une deuxi&me partie d*une deuxi&me 
phase (B), la deuxiime partie 6tant constitute du mat6riau k caract^riser, Ttehantillon 
(L) ayant au moins une dimension (t) suffisamment faible pour permettre un 
rel&chement des contraintes dans ledit tehantillon ; 

10 b. la mesure, sur ledit 6chantillon (L) d*au moins un param^tre de 

deformation (p) d'au moins ladite premiere phase (A), en correspondance d*une 
plurality de points situ6s k des distances difKrentes d*une interface entre lesdites 
premiere (A) et deuxifeme (B) phase ; et 

c. la determination, k partir au moins dudit paramfetre de 

15 deformation (P), d'au moins un parametre mecanique de ladite deuxifeme phase. 

2. Procede selon la revendication 1 , comportant : 

i. la realisation d'une plxuralite d'echantillons (L) qui difBrent 
entre eux pour au moins une propriete geometrique ; 

ii. la mise en oeuvre de retape b. sur chacun desdits echantillons 

20 (L) ; et 

iii. Tutilisation dans retape c. des mesures effectuees sur ladite 
pluralite d'echantillons. 

3. Procede selon la revendication 1 ou 2, dans lequel, pour au 
moins un echantillon (L), retape b. est repetee k au moins deux temperatures 

25 differentes. 

4. Precede selon I'une des revendications precedentes, dans lequel 
retape c* comporte : 

i. la modeiisation du rel&chement des contraintes dans ledit 
echantillon (L) en utilisant une premiere estimation d'au moins une propriete 

30 mecanique du materiau de ladite deuxieme phase (B); 

ii. la comparaison des resultats des mesures de retape b. avec 
ceux de ladite modeiisation ; 

iii. la modification de ladite estimation d'au moins une propriete 
mecanique du materiau de ladite dexjxieme phase et la reiteration des sous-etapes i. k 

35 iii. jusqu'^ la minimisation de I'ecart entre lesdits resultats des mesures et de la 
modeiisation. 
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5. Proc^d^ selon la revendication 4, dans lequel la moddisation est 
une simulation num6rique par ^l^ments finis. 

6. Proc€d6 selon Tune quelconque des revendications pr^c^dentes^ 
dans lequel ledit syst^e composite est choisi pamii : un substrat pr^sentant sur sa 

S surface une couche continue, un substrat prdsentant sur sa surface des ilots ou des 
bandes de metallisation, une couche avec une zone incluse dans le substrat, un 
transistor, une couche & rint^rieur d'un substrat, une matrice contenant des inclusions, 
des fibres ou des filaments. 

7. Proc6d6 selon Tune quelconque des revendications pr^c^dentes, 
10 dans lequel ledit ^chantillon (L) comporte au moins une dimension (t) & T^chelle 

microm^trique ou nanom^trique. 

8. Proc6d6 selon Tune quelconque des revendications pr^c^dentes, 
dans lequel ledit 6chantillon (L) est une lame ayant deux faces sensiblement parall^les 
et dispos^es sensiblement perpendiculairement par rapport k Tinterface entre lesdites 

1 5 premiere (A) et deiixidme (B) phase. 

9. Proc6d6 selon les revendications 2 et 8, comportant la 
realisation d'une plurality de lames (L) d'^paisseurs diff^rentes. 

10. Proc6d6 selon I'une quelconque des revendications 1^7, dans 
lequel ledit ^chantillon est ime lame dispos^e selon un angle par rapport h Tinterface 

20 entre lesdites premiere (A) et deuxieme (B) phase. 

11. Proc6de selon les revendications 2 et 10, comportant la 
realisation d'lrne plurality de lames dispos^es selon des angles diff^rents par rapport k 
rinterface entre lesdites premiere (A) et deuxieme (B) phase. 

12. Procede selon Tune quelconque des revendications pr^cedentes, 
25 dans lequel ledit echantillon est une lame en forme de coin, pr^sentant deux faces 

fomiant un angle entre elles. 

13. Procede selon la revendication 12, comportant la realisation 
d'une pluralite de lames presentant deux faces formant entre elles des angles 
differents. 

30 14. Precede selon Tune quelconque des revendications precedentes, 

dans lequel les mesures prevues k retape b. sont effectuees par diffraction d*\m 

faisceau eiectronique convergent. 

15. Procede selon la revendication 14, dans lequel retape b. 

comporte Tobservation de lignes de Holz pour au moins un plan cristallographique de 
35 ladite premiere phase (A) et la detemiination d'au moins un parametre parmi : la 

largeur des dites lignes de Holz, leur position et leur structure interne. 
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16. Proc^d^ selon la revendication IS, dans lequel T^tape b. 
comporte la determination d'au moins la largeur d*au moins certaines desdites lignes 
de Holz et le calcul, pour chacune d'entre elles, d*une rotation maximale pmax le long 
de I'axe du faisceau d'^lectrons. 

5 17. Proc6d6 selon la revendication 16, dans lequel T^tape c. met en 

oeuvre le trac^ d'au moins une courbe repr^sentant une dite rotation maximale en 
fonction de la distance par rapport k I'interface entre lesdites premiere (A) et 
deuxi^e (B) phase. 

18. Precede selon la revendication 17, dans lequel T^tape c. met 

10 dgalement en oeuvre par simulation le trac6 de courbes repr^sentant la rotation 
maximale Pmax en fonction de la distance par rapport k IMnterface entre lesdites 
premiere (A) et deuxifeme (B) phase pour des valeurs possibles du module d'Yoimg 
et/ou du coefficient de Poisson du mat^riau de ladite deuxifeme phase (B), ainsi que la 
minimisation de T^cart entre les coxirbes simul^es et les courbes exp6rimentales pour 

15 determiner le module d'Yoimg et/ou du coefficient de Poisson du mat^riau de ladite 
deuxi^me phase (B). 
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